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Introduction 

During the past year the City College group continued 
to work in close collaboration with the Goddard Institute 
for Space Studies (GISS) on an evaluation of the feasi- 
bility of long-range weather prediction thi-ough the use 
of global general circulation models (GCMs). The major 
emphasis of the study has been on the quality of simula- 
tions of the- monthly mean state of the atmosphere gener- 
ated by models from real global initial data at the be- 
ginning of the month. Early in the year the last of a 
series of monthly mean simulations was completed v/ith the 
original GISS GCM (Somerville et al., 1974) » after which 
attention was transferred to a new "climate model" de- 
veloped by the atmospheric sciences group at GISS under 
the direction of James Hansen (Hansen et al,, 1978). 

A series of experiments was carried out with an early 
version of the new model, using data for the period Octo- 
ber 1976 through February 1977, to determine the credi- 
bility of global monthly mean simulations generated by a 
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coarse-resolution (8° of latitude by 10° of longitude)^ 
form of the model. The results of these experiments, as 
well as those with the 6ISS GCM, have already been des- 
cribed in a series of technical reports by members of the 
project staff, and only a brief summary is presented in 
this annual report. 

As in the past, the primary motivation of the project 
has been to investigate the feasibility of long-range 
weather prediction with global GCMs. The climate model 
v/as designed for the purposes of very long term climate 
simulation and to test theories of climate change, and 
was not intended for use as an operational v;eather fore- 
casting system. However, as it contains the same dynami- 
cal and physical ingredients as most numerical v/eather 
prediction (N¥P) models and GCMs, it is a suitable vehicle 
for studies of long-range predictability. Furthermore, 
the new coarse-mesh model is very fast. In its most 
recent form, the model generates a one-day global simula- 
tion on the 8 X 10 grid in four minutes (on an IBM 360/95 

^ In the terminology of the GISS modeling group, the 8° x 
10° grid is sometimes identified as "mediiuD”, and the 
terms "fine" and "coarse" may be used to refer to 4° x 5° 
and 12° x 15° grids, respectively. However, in the con- 
text of this report it is appropriate to refer to the 
8 X 10 model as "coarse-mesh" or "coarse-resolution" in 
order to distinguish these experiments from those con- 
ducted with the 4x5 GISS (Somerville et al., 1974) GCM, 
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computer), so that a 30-day forecast can he executed in 
two hours. This is an order of magnitude faster than the 
4x5 GISS (Somerville) model and permits that many more 
experiments to be conducted in comparable time. 

Of course, the high speed of the new model is achieved 
mainly at the price of its coarse resolution, which re- 
quires certain parameterizations of surface boundary 
conditions, as well as inherent filtering of smaller- 
scale features of the initial state, and this may degrade 
the quality of the simulations. Real data tests of the 
model such as the monthly mean prediction experiments 
carried out on this project help to provide realistic 
estimates of model credibility. If a model is to be used 
to investigate the causes of climate change and climatic 
anomalies, it is appropriate to ask how well it can simu- 
late ‘the time-averaged state of the real atmosphere for 
the next month. In this sense, the monthly mean simula- 
tion tests may be viewed as part of the model development 
program. 

Publications during the past year included; 

Spar, Jo, 1977: A summary of monthly mean simulation 

experiments with the GISS model (GSFC). Third NASA 
Weather and Climate Program Science Review, Novem- 
ber 29-30, 1977, NASA, Goddard Space Plight Center, 
Greenbelt, Md, NASA Conference Publication 2029, 
pp. 323-527. Paper Wo. 58, 

Spar, J., J. J, Notario, and W. J. Quirk, 1978: An 

initial state perturbation experiment with the GISS 
model, Mon, Wea, Rev, , 106, 89-100. 
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The following technical reports were distributed: 

Lutz, R. J., 1978: Experiments in monthly mean simula- 

tion of the atmosphere with a coarse-mesh general 
circulation model. The City College, N. Y. , N, Y. , 

46 pages. (M.A. thesis) 

Notario, J. J., 1978: The influence of random initial 

state errors on monthly mean simulations with a 
coarse-resolution atmospheric model. The City College, 
N. Y. , N. Y., 46 pages. (M.A. thesis) 

Klugman, R. , 1978: .The influence of initial conditions 
on monthly mean simulations with a global atmospheric 
model. The City College, N. Y, , N. Y. , 42 pages, 

(M.A. thesis) 

Spar, J, and R. lutz, 1978: Simulations of the monthly 

mean atmosphere for Eebruary 1976 with the GISS model. 
The City College, R. Y. , N. Y. , 29 pages. (Accepted 
for publication in the Monthly Weather Review , ) 

Spar, J,, R. Klugman, R. J. lutz, and J. J, Notario, 

1978: Monthly mean simulation experiments with a 

coarse-mesh global atmospheric model. The City 
College, N. Y. , N. Y. , 56 pages, (Abbreviated version 
submitted for publication to Monthly Weather Review . ) 

A paper on "Monthly Forecasting Experiments with 

Atmospheric Models" was presented in 28 March 1978 by 

J, Spar at the Nev/ York Academy of Sciences, 

Three graduate student research assistants — 

Robert J. Lutz, Jesus J. Notario, and Robert Klugman — 

completed their master's programs at The City College and 

resigned from the project staff this past summer. They 

v/ere replaced in September by tv/o new graduate assistants: 

Zaphiris Christidis and Ronald Filadelfo, Both Notario 

and Klugman have remained at GISS (Lutz left to pursue 

a doctoral program in Maryland), so that there has been 
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no break in the continuity of the project. 


Monthly Mean Simulationg for February 1976 
with the GISS Model 

(a more complete paper on this experiment has 
been distributed as a Technical Report and 
accepted for publication in the Monthly 
Weather Review . ) 

Three monthly mean simulations of the global atmos- 
phere were computed for February 1976 with the GISS model 
(Somerville et al.p 1974) from observed initial conditions 
on the first day of the month. The first two simulations 
(designated TCL and M2) v/ere part of a replication experi- 
ment in which the identical program, with identical ini- 
tial and boundary conditions, was repeated on the IBM 
360/95 computer, the two runs differing only in the 
schedule of interruptions and restarts necessitated by 
other demands on the computer during the month-long fore- 
cast, A comparison of these two forecasts provides a 
measure of the reproducibility of results. Table 1, show- 
ing the results of the replication experiment in terms of 
predicted and observed mean zonal and eddy energies, in- 
dicates that the differences between simulations are 
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Table 1 « Zonal available potential energy zonal 

kinetic energy (Kjyj)» and eddy available 
potential energy (Pg) and eddy kinetic energy 
(Kj,) of standing waves only, for February 1976 
over the Northern Hemisphere and the globe, 
integrated up to the 120 mb level for fore- 


casts Ml and M2 and 
monthly atmosphere. 


Northern Hemisphere 



Ml 

M2 

0 

% 

58.9 

58.0 

49.8 

% 

8.56 

8.42 

7.21 


2.35 

2.45 

3.71 


1.14 

1.51 

2.25 


the observed (0) mean 
Units: 10^ J m“^. 


Globe 

Ml 

M2 

0 

42.5 

41.8 

35.9 

6.69 

6.61 

5.97 

1.69 

1.71 

2.29 

0.91 

1.00 

1.54 
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relatively small compared v/ith the simulation errors. 
Root-mean-square (rms) differences and SI comparison 
scores between the two simulations are shown in Table 2 
for three synoptic fields and for various geographic 
regions. The rms errors of replication over the Northern 
Hemisphere of the monthly mean fields are seen to be ap- 
proximately 2 mb, 20 m, and 1° K for sea-level pressure, 
500 mb height, and 850 mb temperature, respectively, in- 
dicating roughly the computational '‘noise" level of the 
model. 

Table 2 . Root -me an- square (rms) differences and SI 
comparison scores for February 1976 fore- 
casts ML vs. M2. 



Region 

Sea-level 

Pressure 

500 mb 
Height 


850 mb 
Temperature 

rms (mb ) 

SI 

rms (m) 


rms CK; 

1. 

Globe 

1.5 

58 

16 

21 


2. 

Northern Hemisphere 

1.8 

45 

19 

25 

1.2 

5. 

Tropics 

0.8 

57 

5 

55 


4. 

E, Pacific-U. S. 

2.1 


28 


1.2 

5. 

North America 

2.5 

51 

26 

20 


6. 

United States 

2.2 

48 

28 

16 

1,4 

7. 

Europe 

2.1 

45 

19 

25 
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The monthly mean simulation skill for February 1976 
is generally consistent with that found in three earlier 
winter month forecasts with the GISS model, as shown in 
Table 3. Small but consistent skill relative to clima- 
tology is indicated over the Northern Hemisphere for the 
500 mb heights and 850 mb temperatures, but not for the 
sea-level pressure field. 

Table 5 . Summary of rms errors and SI skill scores for 
four GISS-model simulations of monthly mean 
sea-level pressure, 500 mb height, and 850 mb 
temperature over the Northern Hemisphere. 

M denotes the model simulation and C repre- 



sents 

a "forecast” 

of climatology. 






Jan. 

1973 

Jan, 

1974 

Jan. 

1975 

Feb. 

1976 

4 


M 

C 

M 

C 

M 

C 

M 

0 


A. 

rms 

errors 







Sea-level Pressure 

(mb) 

10.0 

8.7 

8,6 

9.2 

5.5 

6.6 

8.8 

6.3 

500 mb Height (m) 


72 

94 

80 

108 

62 

82 

81 

99 

850 mb Temperature 

(K) 

4.1 

4.3 

4.7 

5.1 

4.1 

4.5 

4.4 

4.6 


B. 

SI scores 







Sea-level Pressure 


81 

81 

79 

89 

64 

73 

75 

81 

500 mb Height 


45 

55 

53 

60 

42 

52 

43 

53 
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In a third simulation for February 1976, to measure 
the influence of sea-surface temperature (SST) anomalies 
on the monthly mean forecasts, observed monthly mean SSTs 
were inserted in place of climatological monthly mean 
values as surface boundary conditions. The result was a 
small beneficial impact on the global and hemispheric 
errors, but not on the regional scores. Over the Northern 
Hemisphere, the rms errors for sea-level pressure, 500 mb 
height, and 850 mb temperature were reduced from 8.8 to 
7.2 mb, from 81 to 66 m, and from 4.4 to 4.5° K, respec- 
tively, while the corresponding SI scores for the first 
two fields changed from 75 to 72 and from 45 to 41, re- 
spectively. The effect of the anomalous SST field on the 
simulated synoptic patterns was hardly discernible to the 
eye. ^ 

Monthly Mean Simulation Experiments 
with a Coarse-Mesh Atmospheric Model 

(a more complete paper on this subject has been 
distributed as a Technical Report. The three 
Technical Reports by Lutz, Notario, and Klugman 
also contain additional details on the experi- 
mental results summarised briefly belov/. ) 
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The period October 1976 through February 1977 was 
chosen for a test of an early version of the new coarse- 
resolution (8° x 10°) climate model developed at GISS 
(Hansen et al,, 1978), This was a very anomalous winter 
over North America, with abnormally cold weather in the 
eastern United States and high temperatures in the west. 
Fo\ir groups of experiments were carried out, all based 
on global data for the period provided by the National 
Center for Atmospheric Research (NCAR) and the National 
Meteorological Center (NMG) and derived from operational 
NMC analyses. 

The version of the coarse-mesh model employed in the 
experiments described belov/ (version 252) was derived 
from the GISS model, and employs the same 'Ulatsuno TASU” 
extrapolation scheme as in Somerville et al, (1974), with 
a 12-minute time step. Hov/ever, a number of changes 
have been made in the model physics, notably in the cal- 
culation of radiative transfer, as well as in the param- 
eterization of convection and surface boundary fluxes 
(Hansen et al., 1978). The model is still undergoing 
development, and the experiments with version 252 were 
conducted as part of a program of monitoring progress 
rather than as a final evaluation of the model. 

For the first experiment, the 8° x 10° climate model 
was initialized v/ith global data for 00 GMT on the first 
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day of each of the five months. Climatological monthly 
average SSTs and sea ice fields were used as surface 
boundary conditions. The 12-hourly outputs of the model 
were averaged at the end of each of the five month-long 
forecast runs to produce a set of five monthly mean simu- 
lations. The evaluation of the model simulations was 
limited mainly to the synoptic fields of sea-level pres- 
sxire, 500 mb geopotential height, and 850 mb temperature, 
and expressed numerically in terms of rms errors and SI 
skill scores. To provide a standard for the evaluation 
of the model simulations, monthly climatological fields 
of the three variables, derived from KCAR data, v;ere also 
evaluated as "forecasts" of the five monthly mean states. 

Examination of the observed and simulated monthly 
mean maps (not reproduced here) reveals that the model 
failed to simulate adequately the sea-level pressure 
field. At the 850 mb level, the model exhibited a cold 
bias at all latitudes, especially in higher latitudes, 
with an average error over the Northern Hemisphere of 
- 3.5^ C for the five months. This is reflected hydro- 
statically in low geopotential heights in the simulations 
of the 500 mb level, v/hich average 93 m too low over the 
Northern Hemisphere. The phase opposition betv;een the 
cold east and warm west over North America was also 
poorly simulated by the model, as was the amplitude of 
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the anomalous long v/ave pattern at 500 mb associated with 
the severe North American winter. Rms errors and SI 
skill scores of the model (M) over the Northern Hemisphere, 
together with those for a "forecast" of climatology (C), 
are shov/n in Table 4, where the model simulations are 
seen to be inferior to climatology. 

At this stage of the model’s development, it is ap- 
parently not yet capable of duplicating realistically 
significant departures from climatology of monthly mean 
synoptic patterns. A comparison of the January and Feb- 
ruary Northern Hemisphere results in Table 4 with the 
corresponding error scores for the GISS model, shown in 
Table 3» indicates that the coarse-mesh model simulations 
appear to exhibit even less skill relative to climatology 
than the 4° x 5° GISS model. (However, as the results 
are for different years, they may not be comparable.) 

In a second simulation experiment with the coarse- 
mesh model, the climatological monthly mean sea-siirface 
temperatures (SSTs) were replaced with observed monthly 
mean SSTs derived from satellite radiometer measurements 
(Brower et al., 1976) and provided on a latitude- 
longitude grid from which four point averages were com- 
puted for each coarse-mesh grid point. The magnitude 
(in °C) and distribution of the SST aiiomalies (observed - 
climatology) is illustrated for the month of January 1977 
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Table 4 . Root-mean-square (rms) errors and SI skill 


scores of monthly mean model simulations (M) 
and of climatology (C) over the Northern Hemi- 


sphere for 

winter 

1976-1977 

• 



A. rms errors 

Oct. 

Nov, 

Dec, 

Jan. 

Feb. 

Sea-level pressure (mb) 

M 

5.3 

7.2 

6.1 

7.2 

7.8 

C 

3.0 

2.9 

3.5 

5.9 

4.4 

500 mb height (m) 

M 

104 

114 

112 

131 

119 

C 

39 

40 

43 

73 

55 

850 mb temperature (°K) 

M 

5.1 

5.2 

5.5 

5.9 

5.4 

4- 

c 

2.2 

2.1 

1.9 

2.7 

2.3 

B, SI skill scores 





Sea-level pressure 

M 

95 

102 

90 

89 

87 

C 

65 

53 

51 

64 

62 

500 mb. height 

M 

52 

71 

67 

67 

62 

G 

47 

42 

43 

57 

48 

850 mb temperature 

M 

66 

65 

67 

67 

62 

C 

49 

47 

39 

51 

43 
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in Eig. 1 in the form of a digital global array. (The same 
data source was used for the GISS model SST experiment 
described earlier.) 

The response of the model atmosphere to SST anomalies 
is a complex non-linear combination of local and remote 
effects resulting from the parameterization of the surface 
fluxes of heat and moisture in terms of air-sea tempera- 
ture differences. Locally, we find over positive SST 
anomalies ("hot spots") higher 850 mb temperatures, lower 
sea-level pressures, and slightly higher 500 mb heights 
in, the model simulations with the observed SSTs than in 
the runs with climatological SSTs. However, as shown in 
Table 5 by the rms differences and SI comparison scores 
between the two simulations — A, v/ith the observed SSTs 
and M, with the climatological SSTs — the general impact 
of the SST anomalies on the monthly mean simulations is 
relatively small over the Northern Hemisphere, especially 
at the 500 mb level where it is virtually negligible. 
Furthermore, a comparison of the error scores for the A 
and M simxilations, shown in Table 6 for the Northern 
Hemisphere, reveals no beneficial impact v/hatsoever re- 
sulting from the use of observed SSTs. V/hether this is 
due to poor quality of the SSTs or insensitivity of the 
model, or both, is not known. As noted earlier, the 
GISS model experiment did indicate a small beneficial 
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TAU a .57JOO.O 


ANOMALY SST S FO» JAN 1977 ( CaS . -CL IM . > . DEGHEES KELVIN 


DAY 155A . HOUR 1 2.0 



- 1 no 


- 

150 


- 

120 


-90 



-60 


-30 



0 



30 



60 



90 



1 20 



1 50 





J 

l 

2 

3 

4 

5 

6 

7 

H 

9 

10 

1 1 

12 

13 

1 4 

1 5 

16 

1 7 

IQ 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

LAT 

MEAN 

24 





































90.0 


23 





































82.2 


22 





































74.3 


21 











0 





1 

- 

-1 

0 

0 

















66.5 

0 

20 

0 

0 


-1 

-2 





-1 

2 


3 

1 

-1 

-2 

0 

0 

0 


1 















-1 

53.7 

0 

19 

-3 ■ 

-2 

- 2 

- 1 

- 1 

0 








-1 

6 

-2 

-1 

0 
















0 

-2 

-2 

50.9 

-1 

1 a 

- 4 

-2 

- 1 

- 1 

0 

1 







7 

1 1 

4 

0 

-1 



2 


2 

1 

6 









9 

5 

2 

-6 

43.0 

2 

1 7 ' 

0 . 

0 

0 

1 

1 

0 






6 

2 

1 

1 

2 

1 

2 



1 

1 










4 


3 

2 

0 

35.2 

1 

16 

1 

1 

1 

1 

1 

0 

0 



1 


2 

2 

2 

2 

1 

1 















2 

1 

0 

1 

0 

27,4 

1 

15 

- 1 

0 

0 

0 

0 

- 1 

- 1 

2 



0 


0 

1 

1 

1 

0 






2 


1 

3 


1 



0 

0 

0 

0 

-1 

-1 

19.6 

0 

14 

-3 

- 3 

- 2 

-2 

- 3 

-1 

0 

-1 

0 

1 

-1 


0 . 

- 1 

- 1 

-1 

3 








1 

1 


0 


0 

-1 

-2 

-3 

-3 

-3 

-3 - 

1 1 .7 

-1 

13 

- 2 

- 1 

-2 

— 2 

-2 

-1 

- 1 

-1 

-3 

-1 

-1 



-2 

-2 

-1 

0 

-1 

-2 





1 

0 

0 

1 

0 

0 

0 

0 

-1 

-1 

-1 

-3 

-3 

3.9 

-1 

12 

1 

- 1 

- 1 

0 

0 

0 

0 

0 

0 

-1 






- 1 

-1 

-1 

-2 

-3 




1 

0 

0 

- 1 

-1 

0 

- 1 


0 


-3 

-2 

-1 

-3 .9 

-1 

1 1 

-3 

- 1 

0 

0 

0 

0 

a 

-1 

0 

1 

1 





0 

0 

0 

-1 

0 




0 

-1 

0 

- 1 

-1 

- 1 

0 

1 

2 

1 

-1 

-2 

-2 

- 11.7 

0 

10 

0 

- 1 

- 1 

- 2 

-2 

0 

0 

1 

t 

0 

0 





0 

0 

0 

0 

0 



0 

1 

1 

0 

0 

0 

- 1 

1 




0 

-1 

-1 

- 19.6 

0 

9 

0 

0 

0 

0 

-2 

0 

0 

1 

2 

2 

1 




0 

0 

1 

2 

1 

0 



1 

2 

2 

2 

2 

2 

2 

1 





0 

0 

- 27.4 

1 

8 

0 

1 

1 

2 

I 

-1 

-2 

- 2 

0 

2 

3 



3 

0 

1 

0 

-1 

-1 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

2 

2 

2 



0 

-1 

- 35.2 

0 

7 

2 

1 

1 

1 

0 

- 1 

-2 

0 

0 

-1 

1 


2 

2 

-1 

J 

-1 

1 

1 

-2 

-3 

-1 

-3 

-4 

0 

1 

1 

- 1 

0 

0 

0 

0 

0 

0 

0 


- 43.0 

0 

6 

0 

0 

1 

0 

0 

0 

1 

1 

1 

-1 

-1 


0 

-1 

-2 

-1 

-1 

-1 

-1 

-2 

-2 

0 

0 

- 1 

-2 

- 3 

-3 

-1 

0 

- 1 

- 1 

1 

1 

1 

0 

0 

- 50.9 

0 

5 

1 

0 

0 

- 3 

-2 

- 1 

2 

3 

3 

2 

1 

1 

-1 

-1 

0 

0 

0 

1 

3 

0 

0 

0 

-1 

- 1 

- 1 

• 1 

-2 

-1 

-2 

1 

-1 

1 

1 

-3 

-4 

0 

N 

• 

1 

0 

4 

-2 

-2 

- 1 

0 

- 1 

-1 

-1 

-1 

0 

0 

-1 

-2 


-2 

1 

2 

3 

2 

2 

1 

-1 

-2 

-1 



-1 








-1 

-2 

-2 

- 66.5 

-1 

3 

0 

- 1 

0 

- 1 

-1 









-2 

0 

0 




















0 

- 74.3 

-1 

2 





































- 82.2 


1 





































- 90.0 


J 

I 

' 2 . 

3 

4 

5 

6 

7 

a 

9 

1 0 

1 1 

1 2 

13 

1 4 

1 5 

16 

17 

IB 

19 

20 

21 

22 

23 

24 

25 

26 

27 

26 

29 

30 

^1 

32 

33 

34 

35 

36 

LAT 

MEAN 


- 1 SO 
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150 


- 
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-90 



-6 0 


-30 
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30 



60 



90 



1 20 
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rig. 1 Sea-surface temperature anomaly (°K), January 1977. 


Table 5. 


Root-raean-square (rms) differences and SI com- 
parison scores between monthly mean simulations 
computed with observed (A) versus climatological 
(M) sea-surface temperatures. October 1976 
through February 1977. Northern Hemisphere 
only. 


rms difference 

Sea-level pressure (mb) 
850 mb temperature (°G) 
500 mb height (m) 

SI comparison score 
Sea-level pressure 
850^ mb temperature 
500 mb height 


1976 1977 


Oct. 

Nov . 

Dec. 

Jan. 

Feb. 

1.6 

1.7 

2.4 

2.0 

2.1 

1.6 

1.5 

2.1 

1.8 

1.7 

12 

17 

24 

24 

21 

35 

38 

38 

41 

35 

28 

28 

33 

31 

25 

15 

18 

22 

22 

20 
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Table 6 


. Root-mean- square (rms) errors and SI skill 
scores of monthly mean simulations computed 
with observed sea-surface temperatures (A) com- 
pared with errors (from Table 4) of simulations 
computed with climatological SSTs (M) . October 


1976 through February 1977. Northern Hemisphere, 

1976 1977 





Oct, 

Nov. 

Dec . 

Jan, 

Feb. 

rms errors 








Sea-level pressure 

(mb) 

A 

5.5 

7.4 

6.6 

7.4 

8.2 



M 

5.5 

7.2 

6.1 

7.2 

7.8 

850 mb temperature 

(°0) 

A 

5.2 

5.5 

5.5 

6.0 

5.6 



M 

5.1 

5.2 

5.5 

5.9 

5.4 

500 mb height (m) 


A 

105 

112 

112 

155 

116 

• , 


M 

104 

114 

112 

151 

119 

SI skill scores 








Sea-level pressure 


A 

98 

104 

92 

92 

95 



M 

95 

102 

90 

89 

87 

850 mb temperatiire 


A 

72 

69 

69. 

75 

68 



M 

66 

65 

67 

67 

62 

500 mb height 


A 

55 

71 

70 

76 

61 



M 

52 

71 

67 

67 

62 
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of the observed SSTs. 

A third group of coarse-mesh model experiments was 
carried out to investigate the effects of systematic al- 
terations in initial conditions on the monthly mean simu- 
lations for October 1976 and January 1977. In one set of 
computations, the month-long run was started with data 
for the second and third day of the month, respectively, 
instead of the first, to test the sensitivity of the 
monthly mean simulations to the arbitrary choice of ini- 
tialization time. In still another computation, the 
model v/as initialized with time-averaged initial condi- 
tions based on the first five 12-hourly data sets of each 
month, while in the final experiment of this group the 
model v;as re-initialized for the monthly run v/ith the 
average of the first five days of model output. The 
latter two computations were carried out to test a hypo- 
thesis that smoother initial conditions would produce 
more realistic monthly mean simulations. 

Shifts of one and tv/o days in the initialization time 
were found to produce only small, but not trivial, changes 
in the monthly mean syiaoptic fields. The rms errors and 
SI skill scores of the October 1976 and January 1977 
simulations computed from shifted initial conditions are 
similar to those shov/n in Table 6. Over the Northern 
Hemisphere, the impact of shifting the initialization 
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time, as expressed by the rms differences between simula- 
tions from shifted and unshifted initial conditions, was 
found to be (averaged over the two months) 1.8 mb, 1.5° C, 
and 19 m for a one-day shift, and 2.5 mb, 2.0° C, and 29 m 
for a two-day shift, for sea-level pressure, 850 mb tem- 
perature, and 500 mb height, respectively. These impact 
scores are, of course, much smaller than the simulation 
errors, and may be representative of the inherent minimal 
error of monthly mean simulations associated with the 
arbitrary choice of initialization time. 

The use of time-averaged initial conditions, whether 
of observed fields or of model output, produced no bene- 
ficial effects. The rms errors and SI skill scores were 
found to be essentially unaltered by the use of either 
method of smoothing. The impact of the re-initialization 
method (with time-averaged initial conditions computed 
from the first five days of model output) v/as, in fact, 
negligible. 

The final experiment v/ith the coarse-mesh model was 
a ’’noise-level" estimation for October 1976 based on the 
multiple initial state random perturbation method (Cher- 
vin and Schneider, 1976; Spar et al., 1978), with initial 
global rms errors of 3 ms”^ in wind components, 1° C in 
temperatxzres, and 3 cJb in sea-level pressures. Pour 
perturbations (Pl through P4) of the "control" initial 
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state were generated, yielding a total set of five monthly 
mean simulations. As shown in Table 7» differences among 
the five simulations were virtually negligible over the 
Northern Hemisphere, in terms of rms differences and SI 
comparison scores, compared with the simulation errors 
in Table 6. Average rms differences over the hemisphere 
among the simiilations are approximately 1 mb, 1° C, and 
10 m for sea-level pressure, 850 mb temperature, and 
500 mb height, respectively. These values are smaller 
than those found by Spar ef ai. (1978) in the correspond- 
ing experiment for January 1975 with the 4^ x 5 ^ GISS 
model, which may.be a reflection of either a. less active 
month or the coarser resolution of the climate model. 

. The geographical distribution of the model noise 
level, as represented by global maps of the standard de- 
viations of the simulated variables (not reproduced here), 
shows an increase in noise level from negligible values 
near the Equator to maxim\im values in high latitudes, 
but with considerable zonal variations as well, apparently 
due to different responses over land than over sea. How- 
ever, the influence of 'random initial state errors on the 
monthl 3 ’' mean simulations is insignificant compared with 
the simulation errors, indicating both stable model be- 
havior and the need for further model improvement. 
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Table 7 » Root-mean- square (rms) differences and SI 

comparison scores over the Northern Hemisphere 
between monthly mean simulations for October 
1976, M denotes control simulation. PI, P2, 

P3» and P4 represent simulations from four 
different random perturbations of the control 
initial conditions. 

PI - M PI - P2 PI - P3 PI - P4 

Sea-level rms difference (mb) 0.8 1.0 1.1 1,3 

Pressure SI comparison 23 26 27 29 

850 mb rms difference (°C) 0,9 0.9 0.9 1.0 

Temperature SI comparison 16 18 17 19 

500 mb ‘ rms difference (m) 9.1 9.0 9.4 11.8 

Height SI comparison 10 10 11 14 
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Current Activities and Plans for the Future 


The simulation experiments described above were part 
of a preliminary study of the performance of the climate 
model at an early stage in its development (version 252). 
Since the inception of those experiments, the model has 
imdergone several changes in physics, dynamics, and nu- 
merics (Hansen et al., 1978), At the end of October 1978, 
a new series of monthly mean simulation experiments was 
started for the same five-month period, October 1976 
through February 1977, but v/ith a revised version of the 
model, version 391. As the analysis of results from the 
nev/ computations is still in progress, they are not in- 
cluded in this report, but v/ill be described later. 

Although a description of the climate model lies 
outsi'de the scope of this report, a brief list of changes 
from version 252 to version 591 is given below. 

1. leapfrog extrapolation replaced the TASU-Matsuno 
method (Somerville et al,, 1974), and the time 
step for the 8° x 10° grid was increased from 

12 to 15 minutes. The running time of the 
model has been reduced from eight to four 
minutes per simulated day. 

2. Fourier smoothing is used near the poles. 

3. Changes have been introduced in the schemes for 
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boundary layer fluxes, dry and moist convection, 
radiative transfer, and vertical differencing. 

4. Potential temperature is now a prognostic variable, 

5» Ocean temperatures and ice coverage are inter- 
polated daily from monthly averages. 

6. Vertical diffusion is omitted. 

7. Prediction of snov/, groundv/ater , and ground 
temperature conserves moisture and energy, and 
large-scale precipitation occurs when mean rela- 
tive humidity exceeds 100^. 

In addition to repeating the five monthly mean simu- 
lations with version 391 on the 8° x 10® grid with clima- 
tological surface boundary conditions, we are performing . 
another SST anomaly experiment, but with the observed SST 
valuer processed differently. In the new computation, 
the observed SSTs on the 2-|® grid will be spatially aver- 
aged over the 8® x 10® box represented by each ocean grid- 
point, as was done with the climatological SSI data. Com- 
putations are also being carried out with version 391 on 
a 4® X 5® grid, and the results with the two resolutions 
v/ill be compared. Additional experiments are also planned 
to test the sensitivity of the climate model to variations 
in surface boundary conditions, including snov; cover and 
groiind moisture. The model will also be tested on the 
cold winter of 1977-1978. 
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